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Unidirectional nonreciprocal transport is at the heart of many fundamental problems and applications in 
both science and technology. Here we study the novel design of wave diode devices by engineering 
asymmetric shapes of nonlinear materials to realize the function of non-reciprocal wave propagations. We 
first show analytical results revealing that both nonlinearity and asymmetry are necessary to induce such 
non-reciprocal (asymmetric) wave propagations. Detailed numerical simulations are further performed for 
a more realistic geometric wave diode model with typical asymmetric shape, where good non-reciprocal 
wave diode effect is demonstrated. Finally, we discuss the scalability of geometric wave diodes. The results 
open a flexible way for designing wave diodes efficiently simply through shape engineering of nonlinear 
materials, which may find broad implications in controlling energy, mass and information transports. 

The pursuit of novel devices that are able to manipulate energy, mass and information transmission has 
stimulated huge interests in widespread physical branches such as phononics 1 , photonics 2 " 4 and biophysics 5 . 
As one of the most fundamental and applicable devices, the wave diode can rectify non- reciprocal wave 
propagation, in analogy to electronic p-n junctions, where the transmission coefficient is significantly direction- 
dependent. As such, the device conducts in the forward direction but insulates in the backward one. There are 
many interesting theoretical proposals for thermal diodes 6 " 10 , optical diodes 11 ' 12 , spin Seebeck diodes 1314 and 
acoustic diodes 1516 , and many of them have been successfully verified by experiments 2 " 417 " 21 . 

It has been well established that linear structures are not able to break the reciprocity in reflection-transmission 
once the time-reversal symmetry is preserved 22 . Therefore, both nonlinearity and some kinds of symmetry 
breaking mechanism are essential for the real non-reciprocal propagation. In particular, the optical isolators 
or diodes 23 have been realized by harvesting the properties of nonlinearity and asymmetry 24 " 29 . The enhancement 
of nonreciprocal transmission via cascading has also been studied 30 . And this nonreciprocal transmission has 
even been achieved in PT systems 31 ' 32 . Recently, an interesting design of wave diode is proposed 33 , based on the 
nonreciprocal transmission in a one- dimensional (ID) nonlinear layer structure 1112 . The model system is 
described by the discrete nonlinear Schrodinger (DNLS) equation, which is a reasonable approximation for 
the wave evolution in layered photonic or phononic crystals 34 . The spatially varying coefficients break the mirror 
reflection symmetry and this non-homogeneity generates the asymmetric wave propagation together with non- 
linearity. However, the inhomogeneity of coefficients is not the only way to achieve non-reciprocal wave pro- 
pagation and also it is usually very hard to design the inhomogeneous coefficient in real materials. 

To circumvent the difficulty, in this article we report another novel design by seeking the symmetry breaking 
from pure geometric point of view, i.e, engineering asymmetric shapes of homogenous materials. The advantage 
of utilizing geometric asymmetry is guaranteed by the easy and convenient fabrication of functional devices 
merely through shape engineering. Different functionality can thus be tunable by simply tailored adjusting of the 
geometric shape of the target device. Therefore, the design of a wave diode device purely induced by geometric 
effect is very important in the practical point of view and would have found potential applications in the future. 

To unraveling the underlying mechanism of the nonreciprocal geometric wave diode, we first analytically solve 
the transmission problem for a simple model, where the functional nonlinear part consists of a triangle structure 
coupled to a linear ladder lattice in the left and a linear ID chain in the right. The derived transmission coefficients 
exhibit clear asymmetry depending on the direction of the incident wave, thanks to the nonlinear triangle 
structure. The strength and sign of the non-reciprocity can also be tuned by the intensity of the incident wave. 
We then investigate a more realistic system with typical asymmetric geometry. The system is complicated and 
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analytical result is no longer available so that we resort to intensive 
numerical simulations of the propagation of wave packets. The 
results clearly demonstrate non-reciprocal wave propagations 
induced by the asymmetric geometry of nonlinear materials. We 
further verify that our proposed wave diode devices will not change 
the frequency of the incident wave packet which would be important 
in practical applications. Finally, we discuss the scalability of the 
geometric wave diode. 

Results 

Analytic results for a simple geometric diode. Figure 1 (a) shows the 
sketch of a simple geometric wave diode model with non-reciprocal 
wave transmission induced by asymmetric shape. The central 
triangle lattice with L c = 2 interfacial layers is made of nonlinear 
materials and plays the role of functionality for the device. The 
nonlinear triangle is then coupled with semi-infinite linear ladder 
lattice in the left and ID linear chain in the right, respectively. 

The evolution of wave amplitude at each site / follows the time- 
dependent DNLS equation 35 : 

'mi 

mi 



i^ l = {U l + n l )xli l +g l \xli l \ 2 xli l -Y J ^n 



(1) 



where \J\ and n\ depict the on-site energy and the number of nearest- 
neighbors of site /, gi is the nonlinear strength and m/ denotes the 
index of nearest-neighbors of site /. The on-site energy and nonlinear 
strength have homogeneous non-zero values Ui = Uandgi = g only 
in the central functional part. Therefore the symmetry breaking only 
comes from the geometry. For the right ID linear chain, the disper- 
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Figure 1 | (a) Sketch of the geometric wave diode. The central L c = 2 
interfacial layers are nonlinear, with asymmetric geometry. Its left is 
coupled with a semi-infinite linear ladder lattice. The right is connected to 
a semi-infinite ID linear lattice. The linear structures are free of on-site 
potentials and the central 2 layer interfaces are with identical nonlinear 
strength g and on-site potential U. The arrows schematically denote the 
forward incident direction, (b) Transmission coefficients and rectification 
factor as a function of the intensity of transmitted waves. The wave vectors 
are identical with k = 1.57 for both left and right incoming waves. The on- 
site energy U = —3.5 and nonlinear strength g = 1 are set for the central 
two layers. 



sion relation can be obtained as co = 2 — 2 cos k R , where the plane 
wave solution has the form of ^ / = A/e~^ _w) . While for the left 
linear ladder, there are two branches for the energy band. Depending 
on whether bilateral amplitudes on the ladder are (anti-) symmetric 
\jj v = ±xjj l where /' denotes the site in the lower chain of the left 
ladder, the dispersion relations of the plane wave solutions are 
w = 2 — 2 cos k L and co = 4 — 2 cos k L , respectively. 

We can analytically solve the transmission problem of the wave 
propagation by the scattering approach similar as in Ref. [33] . For the 
incident wave coming from left to right, the solutions of the wave 
amplitudes in the central four layers [see Fig. 1(a)] can be expressed 
with the following boundary conditions: 

^, Q = ^=I + R 9 xj/ l = xj/y = Ie ikL + Re~ ikL 

ijj 2 = Te ik \ ^3 = ^2*** 



(2) 



where I, R and T denote incident, reflected and transmitted wave 
amplitudes, respectively. Without loss of generality, we only consider 
the lower dispersion branch of the left ladder lattice, which is ident- 
ical to that of the right ID chain by assuming \jj y = \\i x . As such, both 
left and right linear structures have the same dispersion and the wave 
vectors are identical with k LyR = k under the same energy. Therefore, 
the forward transmission coefficient Gf (where subscript / denotes 
forward) from left to right can be derived as the function of trans- 
mitted wave intensity | T\ 2 (see Methods for the detailed derivation): 



8 sin k L sin k R 



\T\ sink R = 

2|J| 2 sinA; I |2- (a-^)(£-^)| 2 



(3) 



where ft = U + 3 - co + g\T\ 2 andot = U + 2 - co + g\T\ 2 (P 2 + 1 
2/5 cos k R )/4. The ratio sin k R /sin k L appears in the definition as a 
renormalization factor because the group velocities of incident waves 
of the left and right leads might not be the same (if the upper energy 
branch of the left side is chosen, in our case). The factor 2 in the 
denominator is the consequence that there are two incident channels 
from the left due to the ladder structure. One can readily verify the 
conservation condition 2\I\ 2 sin k L = 2\R\ 2 sin k L + \T\ 2 sin k R . In the 
same manner, the backward transmission coefficient G b (where sub- 
script b denotes backward) from right to left can also be obtained as: 



G b :~- 



8 sin k L sin k R 



2\T b \smk L 

\I b \ 2 smk R ~ \2-(a b -e ik ^ b -e^)\ 2 



(4) 



with ot b = U + 2 - w + g\T b \ 2 and P b = U + 3-co + 
g\ T b | 2 (ol\ + 1 — 2d b cos k R ) . One can also readily verify the conser- 
vation in the backward direction that \I b \ 2 sin k L = \R b \ 2 sin k L + 
2 1 T b | 2 sin k R . The factor 2 is due to the fact that there are two trans- 
mitted channels to the left. 

Although the forward and backward transmission coefficients Gf 
and G b have similar expressions, the non-identical behavior between 
parameters (a, /?) and (ot b , p b ) implies the existence of direction- 
dependent transmission coefficients. For plane waves with same 
intensities, the incident waves will be scattered differently by the 
central nonlinear triangle due to its asymmetric shape in respect to 
the forward and backward waves, which leads to the non-reciprocal 
propagation. To quantify the magnitude of the non-reciprocity, a 
rectification factor F R can be defined as: 

Gf — G b 



Fr = 



G f + G b 



(5) 



where a non-zero value of F R emerges for the non-reciprocal wave 
propagation and reaches to its maximal asymmetry with values F R = 
±1. 

We plot Gf and G b as well as F R in Fig. 1(b), as a function of the 
transmitted wave intensities |T| 2 = 2|T^| 2 according to Eq. (3) and 
(4). (Note here, we choose the same transmitted intensities for the 
convenience, because the transmission is analytically expressed in 
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terms of transmitted intensities. We have checked at the stable 
regime that the same incident amplitudes also produce the nonreci- 
procal transmission. In fact, in the following simulations, we choose 
transmissions at the same incident amplitudes for the comparison.) 
We can see that the dependence of G^and G b on transmitted intens- 
ities are almost opposite to each other. As a result, significant non- 
reciprocal effect of wave propagation emerges and the rectification 
factor F R can even change its sign as the modulation of transmitted 
wave intensities. We would like to remark that analytically we obtain 
the nonreciprocal transmission because a 7^ ot b , 7^ fi b that are 
consequences of both the nonlinearity and the asymmetric geometry. 
If the shape is symmetric, we of course have symmetric transmis- 
sions. If the nonlinearity vanishes g — » 0, even with keeping the 
asymmetric geometry, we can see analytically that a = oc b and ft = 
$ b so that the reciprocal transmissions are also recovered. 

Simulation results for a more realistic geometric diode. We now 

consider a more realistic model with typical asymmetric geometry 
displayed in Fig. 2(a). The central functional part with layer index 
1 < / < L c is nonlinear in nature and has asymmetric shape. The same 
left and right linear ladder lattices are coupled with the central part. 
The system is too complicated to be tackled analytically. Therefore, 
numerical simulations will be performed to investigate the wave 
propagation in such a setup. 

We consider the propagation of a pulse or say a Gaussian wave 
packet with the initial condition: 



(6) 



where d and Z 0 is the width and central position of the wave packet, 
respectively. The numerical transmission coefficients Gf and G b are 
defined as the sum of transmitted density norms ^^ 1>L |*A/M| 2 at 
the end of time evolution normalized by the sum of initial density 
norms ^ i \i// l (t = 0)\ 2 .ln Fig. 2(b), we plot the numerical forward 
and backward transmission coefficients Gf and G b as well as the 
numerical F R as a function of incident wave intensities \I\ 2 . 
Although the profiles of forward and backward transmission coeffi- 
cients exhibit similar behaviors, a lag between them is enough to give 
rise to a profound non-reciprocal wave diode effect. 

Figure 3(a,b) displays typical examples of time evolution for the 
same initial wave packets with opposite incident directions. In the 
forward direction, most of the incident waves have been reflected by 
the nonlinear central part. While in the backward direction, most of 
the incident waves can pass through the nonlinear central part. 
Accordingly, figure 3(c,d) shows the distinct density distributions 
within the nonlinear triangular part for forward and backward wave 
transmissions, respectively. The snapshots evidence that the forward 
wave is blocked at the smaller end of the triangle when the wave is 
incident from the left to the right, while the backward wave success- 
fully passes through the central part when incident from the larger 
end of the triangle (right) to the smaller one (left). The interplay 
between nonlinearity and asymmetric geometry within the central 
functional part yields the non-reciprocal wave propagations. The 
asymmetric behavior can be modulated by the incident wave intens- 
ities and wave vector k. More importantly, the sign and magnitude of 
this non-reciprocity can also be manipulated by the geometric shape 
of the central part. Further, to check whether the frequency of the 
wave packet can be shifted or not after the scattering, we plot the 
power spectra of the wave packet at the end of time evolution in 
Fig. 3(e,f). Although nonlinear resonances can be generated by the 
central part, the reflected and transmitted wave packets almost main- 
tain their original frequencies. This effect is important for many 
practical applications since usually we do want to retain the fre- 
quency of the original signal after it passes through the functional 
devices. 
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Figure 2 | (a) Sketch of a more realistic geometric wave diode. The central 
L c layers are nonlinear, with asymmetric shape. The left and right are 
coupled with the same semi-infinite linear ladder lattices. The linear 
structures are free of on-site potentials and the central layers are with 
homogeneous nonlinear strength g and on-site potential U. The arrows 
schematically denote the forward incident direction, (b) Transmission 
coefficients and rectification factor as a function of the intensity of incident 
waves. The initial wave intensities are the same for both directions with 
III 2 = I Iy I 2 . The absolute values of wave vectors are identical with \k\ = 
1.57 for both left and right incident waves. The central part has L c = 6 layers 
with homogeneous on-site energy U = 1 and nonlinear strength g = 1. 
Both of the left and right linear ladders have L/ = L r = 500 layers. The width 
of the initial wave packet is d = 100 and the central position of the initial 
wave packet is l 0 = 250 (750) for left (right) incident wave. The vertical line 
denotes the initial wave intensities III 2 used in Fig. 3. 

Finally, we discuss the scalability of the rectifying function of geo- 
metric wave diodes by linking them in series. As shown in Fig. 4, 
when scaling up to multiple diodes, the rectification is enhanced. This 
can be understood by the following very rough but intuitive reason- 
ing: Assuming n geometric diodes in series are independent and not 
interact with each other, then their forward transmission coefficient is 
n times multiplication of the singe one, expressed as G". Accordingly, 
their backward transmission coefficient will be G b . As such, the 
total rectification factor for the n geometric diodes expresses as 

1 - (^j ^ I (\ + (^j ^ . Readily, it can be proved that 

increasing n will always enhance F^ to its extremum values ± 1. 
However, this is just a rough estimation for ideal cases. In reality, 
the situation is more complicated: the enhancement of the rectifica- 
tion will saturate at some optimal number (5 in Fig. 4), or decrease, or 
even reverse the sign, depending on the wave vector k and intensities 
(not shown). More importantly, scaling up the diode devices will 
steadily suppress the transmission since its value is always smaller 
than 1, i.e., G n decreases with n since G < 1 generally. Therefore, in 
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Figure 3 | (a,b) Time evolution of left and right incoming Gaussian wave packets demonstrates the device as a wave diode. The parameter setup is 
the same as in Fig. 2 with the initial wave intensities III 2 = 0.16 as marked by the vertical line there. (c,d) The density distributions within the 
triangular part for forward and backward cases, respectively. The snapshots are taken at the time when the center of the wave packet is traveling to the 
inside of the triangular part. (e,f) Power spectra of the real part of at the ending time of the evolution for left and right incoming wave packets, 
respectively. The initial wave vectors are I k\ = 1.57, represented by the vertical dotted lines in left and right panels. The ending time of the evolution is 
t = 253. 



practice, one need find the optimal number of geometric wave diodes 
to balance demand of both the transmission and the rectification. 

Discussion 

Geometric phonon and electron diodes have been studied both the- 
oretically and experimentally 38 43 . However, to our best knowledge so 
far there is no existing exact results to clearly expose the underlying 
physics, especially for the more general geometric "wave" diode, 
before this work. In particular, the role of nonlinearity (many-body 
interaction beyond quadric order) to the geometric diode has not 
been appreciated. 

Therefore, it is worthwhile to emphasize here that the asymmetry 
alone in linear (quadric) systems cannot induce non-reciprocal wave 
propagation. Asymmetry must come to play together with nonli- 
nearity to give rise to non-reciprocal wave diodes. In all our studied 
systems, we have either analytically or numerically checked that the 



effect of non-reciprocal wave propagation will totally disappear if the 
nonlinearity is tuned to be zero while keeping all other asymmetric 
settings. These would explain why the electric and thermal rectifica- 
tions measured in the asymmetric graphene and its oxide systems are 
too low 40,42 [both are below 0.13 according to our Eq. (5)], because in 
these graphene- related systems electrons and phonons are weakly 
interacting (especially at low temperatures) so that the nonlinearity 
(high order many-body interaction) is nearly absent. To observe the 
significant rectification in those systems, one need go to the high 
temperature regime, where the nonlinear effect due to many-body 
interactions can be enhanced. (Note that extremely low temperature 
can also produce and enhance the nonlinear effect due to many-body 
interactions, e.g., the nonlinear term in the Gross -Pitaevskii equation 
for Bose-Einstein Condensates.) 

The interplay between some kind of symmetry breaking mech- 
anism and nonlinearity is essential to realizing non-reciprocal wave 
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Figure 4 | (a) Sketch of a geometric wave diode linked in series. The n 
diodes (n = 3 in the sketch) are linked in series to enhance the rectification 
factor. The arrows schematically denote the forward incident direction, (b) 
Rectification factor as a function of the number of diodes linked in series. 
The central part for each diode has L c = 3 nonlinear layers with 
homogeneous on-site energy U = — 1 and nonlinear strength g= 0.1. The 
linear ladder between neighboring diodes has L# = 50 layers and the most 
left and right linear ladders has L\ = L r = 1000 layers. The wave vectors and 
width of initial wave packets are \ k\ = 0.8 and d = 30, respectively. The 
central positions of the initial wave packets are in a distance of 250 layers to 
the nearest end and the ending time is determined when the travelling 
distance of the central peak equals to L — 500 where L is the total number of 
layers of the whole system. 

propagation. Ref. [33] breaks the symmetry by setting spatially vary- 
ing inhomogenous coefficients to induce asymmetric propagation. 
Here we proposed a novel design of the so called geometric wave 
diode by tailoring the asymmetric shape of the central functional part 
of the device. The functional nonlinear part of asymmetric shape can 
be conveniently fabricated in practice and easily be replaced by other 
shapes with different or even reversed rectifications. Therefore, the 
most manipulability of the non-reciprocal wave propagation can be 
solely achieved by geometry engineering of the central functional 
part. 

In conclusion, we have shown the non-reciprocal geometric wave 
diode by engineering asymmetric shapes of nonlinear materials. For 
a simple model, we have obtained analytical results, which explicitly 
prove that the non-reciprocal wave propagation can be induced by 
the asymmetric shape in nonlinear materials. We then have per- 
formed numerical simulations to study a more realistic diode model 
with typical geometric asymmetry. Profound non-reciprocal propa- 
gations of wave packets have been verified with tailored design of 
interplay between nonlinearity and geometric asymmetry. Most 
importantly, the sign and magnitude of the non- reciprocity of the 
wave propagation can be manipulated by shaping the central func- 
tional part into different geometric structures. 

Since the DNLS has found widespread implications in optics, cold 
atoms and spin electronics, the proposed design in this article might 



also find its applications in controlling the light propagation in non- 
linear optical wave guides, cold atom transport in atomtronics 44 ' 45 , or 
spin wave transport in spin caloritronics 14 . In the future, it would be 
interesting to explore how the asymmetric geometry optimizes the 
rectification of energy, mass or information transmission, which may 
have a deep connection to the quantum chaos 46 with nonlinearity. 

Methods 

In order to analytically solve the transmission problem for the first model, we find the 
traveling wave solution of ij/i ~ e~ kot for Eq. (1) as: 



w\j/ l = {Ui + ni)\j/ l +gi\\jf l \ 2 \j/ l -^2\j/ mi 

mi 

Apply the above equation to site 2 in Model I, we obtain: 

e# 2 = (l/ + 3)^ 2 +g\^ 2 \ V 2 - OAi + <Ar + W 



(7) 



(8) 



it can be regrouped by knowing that i/^ = xj/r and = e lkR \j/ 2 from the boundary 
conditions of Eq. (2): 



(j»-«*')*2 



with 



P=U+3-ca+g\T\ 2 
Apply Eq. (7) to site 1, we obtain: 

a4 x = (17 + 3)^ +g\\lf 1 1 Vi - OAo + <Ar + <W 
After regrouping and substitution we can obtain: 

*0=Q«(j&-^)-l)*2 

with 

a = U + 2-co+g\T\ 2 (p 2 + l-2Pcosk R )/4 



(9) 



(10) 



(11) 



(12) 



(13) 



From the first two boundary conditions in Eq. (2), the incident amplitude I can be 
expressed as: 



(14) 



Substitute the expression of \J/ 0 and \J/ l into above formula and noticing that \J/ 2 = Te lkR 
in Eq. (2), we can obtain the relation between incident and transmitted amplitude: 



1 = 



t e- ik L(oL(p-e ik «)/2-l)-(p-e ik *)/2 r 



(15) 



The transmission coefficient of Eq. (3) can thus be derived analytically. For the 
reversed direction, the transmission coefficient of Eq. (4) can also be derived with the 
same consideration. 

For the numerical simulation, we apply the symplectic PQ method introduced in 
the appendix of Ref. [36] to integrate the time dependent DNLS equations. It is a 
SBAB 2 based symplectic algorithm 37 . The dimensionless time step At = 0.02 has been 
used in all the numerical simulations. 
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